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Abstract. We perform a forecast analysis on how well a Euclid-like photometric galaxy 
cluster survey will constrain the total neutrino mass and effective number of neutrino species. 
We base our analysis on the Monte Carlo Markov Chains technique by combining informa- 
tion from cluster number counts and cluster power spectrum. We find that combining cluster 
data with Cosmic Microwave Background (CMB) measurements from Planck improves by 
more than an order of magnitude the constraint on neutrino masses compared to each probe 
used independently. For the ACDM+m^ model the 2a upper limit on total neutrino mass 
shifts from m u < 0.35eV using cluster data alone to m u < 0.031eV when combined with 
Planck data. When a non-standard scenario with iV e g ^ 3.046 number of neutrino species is 
considered, we estimate an upper limit of N e s < 3.14 (95%CL), while the bounds on neutrino 
mass are relaxed to ^ m u < 0.040eV. This accuracy would be sufficient for a 2a detection of 
neutrino mass even in the minimal normal hierarchy scenario C^2m u ~ 0.05 eV). In addition 
to the extended ACDM+m^+iVefr model we also consider scenarios with a constant dark en- 
ergy equation of state and a non-vanishing curvature. When these models are considered the 
error on ^ m v is only slightly affected, while there is a larger impact of the order of ~ 15% 
and ~ 20% respectively on the 2a error bar of iV e ff with respect to the standard case. To assess 
the effect of an uncertain knowledge of the relation between cluster mass and optical richness, 
we also treat the ACDM+m^+A^g case with free nuisance parameters, which parameterize 
the uncertainties on the cluster mass determination. Adopting the over-conservative assump- 
tion of no prior knowledge on the nuisance parameter the loss of information from cluster 
number counts leads to a large degradation of neutrino constraints. In particular, the upper 
bounds for ^ m v are relaxed by a factor larger than two, Yl m v < 0.083 eV (95%CL), hence 
compromising the possibility of detecting the total neutrino mass with good significance. We 
thus confirm the potential that a large optical/near-IR cluster survey, like that to be carried 
out by Euclid, could have in constraining neutrino properties, and we stress the importance 
of a robust measurement of masses, e.g. from weak lensing within the Euclid survey, in order 
to full exploit the cosmological information carried by such survey. 
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1 Introduction 

Over decades neutrino oscillation experiments have provided conclusive evidence that neutri- 
nos have non-zero masses. Such experiments provide constraints on the neutrino mass squared 
difference, while they are not sensitive to the absolute scale of neutrino masses. The latest 
measurements, using solar, atmospheric, and reactor neutrinos, give mass difference between 
neutrino species of Am 2 2 = 7.5 x 10~ 5 eV 2 and lAmf^l = 2.3 x 10 _3 eV 2 [e.g. 1, 2], which in 
turn translates to a lower bound on the sum of the three masses, ^2m u , at 0. 056(0. 095)eV 
in the normal (inverted) hierarchy. On the other hand, cosmological data provides a tool to 
constrain neutrino masses due to the effects neutrinos induce on background evolution and 
growth of structures: relativistic neutrinos affect the Cosmic Microwave Background (CMB) 
anisotropics, whereas, at low redshift when they become non-relativistic, neutrinos suppress 
matter density fluctuations at small scales [e.g. 3, 4, for a review]. 

Given these multiple effects massive neutrinos leave an imprint on many cosmological 
observables; indeed, current neutrino constraints from cosmology rely on a combination of 
data from CMB experiments, Baryonic Acoustic Oscillations (BAOs) measurements, Super- 
novae distance moduli, galaxy clustering and galaxy cluster mass function. Recent constraints 
on the upper limit of the total neutrino mass lie in the range Yl m v < 0.3 — 0.8 eV (95% 
Confidence Level (CL)) [e.g. 5-10], with notable exception of Lyman— a data, which presents 
an even tighter bound of 0.17 eV [11]. 

The number of active neutrinos is known to be three to high precision through the 
measurement of the invisible width of the Z boson at LEP [jV a = 2.9840 ± 0.0082; 12], 
however, the possibility remains that additional "sterile" species exist [e.g. 13]. In fact, in 
order to explain the results of short baseline neutrino oscillations experiments LSDN [14] and 
MiniBooNE [15], as well as the recently discovered reactor neutrino anomaly [16, 17], models 
with one or two light sterile neutrino have been proposed [e.g. 18-20]. Observations of the 
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CMB also seem to point to the same direction, favoring the presence of extra relativistic 
degree of freedom at the time of decoupling, in terms of the effective number of neutrino 
species N e g. For instance, measurements of CMB anisotropies from the South Pole Telescope 
(SPT) [21] and the Atacama Cosmology Telescope (ACT) [22], combined with seven-year 
Wilkinson Microwave Anisotropy Probe (WMAP) data [23] and BAO and Hq measurements, 
have measured iV e ff = 3.71 ± 0.35 and N e g = 3.52 ± 0.39 at 68% CL, respectively, thus 
suggesting values higher than those expected in the canonical scenario (iV e fj = 3.046). 

Among the different probes of the Large Scale Structure (LSS), many works have already 
been proven the ability of galaxy clusters in constraining cosmological parameter, through 
both the evolution of their mass function [e.g. 24-31], and their large-scale clustering prop- 
erties [e.g. 32-34]. Indeed, galaxy clusters supply cosmological information in different ways. 
The evolution of the galaxy cluster population depends on cosmological parameter both 
through the linear growth rate of density perturbations and the redshift dependence of the 
volume element. Furthermore, also clustering of galaxy clusters is sensitive to cosmologi- 
cal parameters through the grow rate of perturbations. These properties makes clusters an 
excellent probe of the growth of structure, and thus a valuable tool to constrain neutrino 
properties [e.g. 29, 35]. 

Recent constraints from galaxy cluster samples combined with CMB data, observation of 
BAO and supernovae Type la are 0.52 eV, N cS = 3.6t\i I 36 ] and 0.72 eV, N cS < 4.6 

(95% CL) [29, 37]. 

While current constraints based on galaxy cluster data rely on relatively small sam- 
ples of clusters identified at redshift below one [e.g. 38, 39], next generation of X-ray (e.g., 
eROSITA 1 , WFXT 2 ), Sunyaev-Zeldovich (e.g., CCAT 3 , SPT-3G) and optical (e.g., DES 4 , 
LSST 5 , PanSTARRS 6 , Euclid 7 ) surveys are expected to increase by orders of magnitude 
the number of galaxy clusters detected, further extending the probed range of redshift up to 
z ~ 2. Such large cluster surveys will provide tight constraints on cosmological parameters, 
independently and complementary to those recovered from other cosmological probes. In 
this work we explore the cosmological information contained in the cluster catalog that will 
be provided by the photometric redshift survey of ESA's Euclid mission, which has been ap- 
proved for lunch in 2019. Specifically, we will make use of cluster number counts and cluster 
power spectrum to derive forecast errors on the total neutrino mass and effective number 
of neutrino species, for a Euclid-like galaxy cluster survey. A similar analysis based on the 
Fisher Matrix approach has been proposed by [40]. Even though the Fisher Matrix technique 
has the advantage of allowing for a quick, analytic estimate of the confidence limits, on the 
other hand it approximates the likelihood function as a multivariate Gaussian function of 
the model parameters. In general, this turns out to be a coarse approximation since the 
likelihood function can be highly non-Gaussian [e.g. 41]; moreover the results obtained with 
this technique depend on the step chosen in the calculation of numerical derivatives with 
respect to the parameters [e.g. 42]. For these reasons we choose to use a more robust forecast 
method based on the sampling of the full likelihood function using a Monte Carlo Markov 
Chain (MCMC) approach. Finally, unlike [40], in our work we explore also the non-standard 

1 http://www.mpe.mpg.de/eROSITA 
2 http://www.wfxt.eu/ 

3 http://www. ccatobservatory.org/index.cfm 
4 http: / /www.darkenergysurvey.org/ 
5 http://www. lsst.org/lsst/ 
6 http : / /pan-st arrs . ifa.hawaii.edu /public/ 
7 http://www. euclid-ec.org/ 
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Figure 1: Effects of the variation of ^ m u and N e g on the linear matter power spectrum 
(left) and halo mass function (right) at z = 0; all other parameters (O m , f^A, Hq, n s , Aj^, r) 
are kept fixed to the WMAP 9-yr best-fit values for ACDM. See text for comments. 



scenario with iV e fj > 3 and ^ m v > 0, to assess the effects of the correlation between the two 
parameters [see e.g. 37, 43]. 

The paper is organized as follows. In Section 2, we briefly explain the physical and 
observable effects of iV e g and m v ■ Section 3 describes the formalism which we use to 
compute the cluster number counts (§3.1) and power spectrum (§3.2), while in § 3.3 we 
outline our forecasting procedure and specify the characteristics of the galaxy cluster survey 
analyzed. Our results for different cosmological model are presented in section 4, and finally 
in section 5 we draw our conclusions. 



2 Massive neutrino effects 



In this section we briefly review the physical effects of N e g and m v an d how they affect 
the halo mass function and the matter power spectrum, two quantities strictly related to the 
observables used to derive the forecast error (see next section). 

What is actually constrained by cosmological data is not the neutrino mass, but the 
neutrino energy density p v , which can be related to the neutrino mass through the relation: 

u ~ Pc ~ 93.14eV/i 2 ' 1 ' 

where p c is the critical density of the Universe, and N v the number of massive neutrinos. 
In the present work, we assume three degenerate massive neutrino states, so that the total 
neutrino mass can be written as ^ m u = 3m u . A degenerate mass spectrum is justified on one 
side by the smallness of the mass difference measured, on the other by the incapability of the 
Euclid cluster data to probe directly the neutrino mass hierarchy, as we will explain later in 
§ 4.2. The effective number of neutrino species parameterizes the contribution of neutrinos 
(or any other relativistic species) to the radiation content in the radiation dominated era 
through the formula: 

7/4 \ 4/3 



Pi 



1 n J NeS 



Pl , (2.2) 
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where p T is the energy density in relativistic species and p 7 is the energy density of photons. 
N e g can be expressed as the sum of the number of massive neutrinos N v and the contribu- 
tion from extra relativistic degrees of freedom, AiV e g. The standard value for the number of 
massive neutrinos is N v = 3.046, where the 0.046 accounts for a non-instantaneous neutrino 
decoupling process and flavour neutrino oscillations effect [44]. Instead, any significant devi- 
ation from AN e g = could indicate the presence of new physics beyond the standard model; 
adding an extra (thermalized) light fermion would contribute AiV e g = 1, but more generally 
a non-integer AiV e ff value could arise from different physics, such as lepton asymmetries [45], 
partial thermalisation of new fermions [19], particle decay [46], non-thermal production of 
dark matter [47], gravity waves [48] or early dark energy [49]. 

Modification of these parameters has effects both on the background evolution and on 
structure formation. Neutrinos decouple from the primeval plasma at T ~ 2 — 4MeV, when 
they are still ultra-relativistic, and given their small mass, ^ m v < leV, they become non- 
relativistic (only) after recombination. During this epoch, their energy density contributes as 
radiation rather than matter, thus changing the expansion rate (via the Friedmann equation 
H 2 = ^p(/9 m + Pr)) and the time of matter-radiation equality (a cq = Pr/(/°m ~~ Pu))- The 
epoch of equality set the time at which sub-horizon density fluctuations start to collapse 
under the action of gravity and structures can evolve. When the other parameters are kept 
fixed, a larger value of ^2 m v or N e s corresponds to a larger value of the radiation density, 
and therefore a postponed time of equality. Such modification is seen in the matter power 
spectrum as a shift to larger scale of its peak, which is determined by the size of the particle 
horizon at the time of matter-radiation equality. Moreover, since on sub-Hubble scales density 
fluctuations grow more efficiently during matter dominated epoch (i.e. after equality), the 
matter power spectrum is suppressed on small scales relatively to large scales (left panel of 
Fig. 1). These effects determine also a suppression of the halo mass function (right panel of 



After thermal decoupling neutrinos constitute a collisionless fluid, whose constituent 
fluid elements free-stream with a velocity, on average, equal to their thermal velocity t>th- 
As the Universe expands, i>th decays adiabatically till neutrinos become non-relativistic. At 
this stage neutrinos behave as hot dark matter particles, suppressing density fluctuations on 
scale smaller than their free-streaming length: 



Due to their high velocity neutrinos cannot be confined on region smaller than their 
free-streaming length, thereby suppressing the density perturbations by a factor proportional 
to (1 — Q u /£l m ). Furthermore, the absence of gravitational back-reaction effects from free- 
streaming neutrinos slow down the grow rate of CDM/baryon perturbations at late times [see 
e.g. 4]. Both these effects are observable in the matter power spectrum as a suppression on 
small scales (left panel of Fig. 1), with a constant amplitude AP(k)/ P(k) ~ —8(Q U /Q m ) at 
k ~ 1/i/Mpc, for linear structure formation [50]. Similarly, the halo number density decreases 
due to the free-streaming neutrino action (right panel of Fig. 1). Again the suppression is 
larger for larger neutrino mass, and for a given ^ m u it is more pronounced for the heaviest, 
late forming haloes, since the damping in power from massive neutrinos shifts the maximum 
cluster mass down (i.e. the scale beyond which the halo mass function is exponentially 
suppressed). 



Fig. 1). 




(2.3) 
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To actually detect the signature of massive neutrinos in the matter power spectrum 
and halo mass function we need a cosmic tracer of such quantities. As mentioned in the 
introduction we will use the cluster number counts and cluster power spectrum as observables 
to detect massive neutrino signature. In the next section we will describe the formalism we 
use to quantify such an effect on the two observable and outline the forecasting procedure 
adopted. 

3 Methodology 

We focus on the galaxy cluster sample which will be provided by a Euclid-like photometric 
survey. For this observable we quantify its clustering statistics in terms of the average 
cluster power spectrum, and its mass distribution in terms of the cluster number counts. In 
describing these two quantities we use the same notation as in [51]. In order to extract the 
cosmological parameter errors from the mock Euclid data, we perform a Bayesian likelihood 
analysis. 

3.1 Cluster number counts 

The number of clusters expected for a survey having sky coverage AO with an observed mass 
between M,* and M^ +1 and redshift between z\ and z\+\ can be expressed as: 

/■ 2 i+i dv r M ° b m+i dM ob r°° k 

Ni, m = AnJ dz — J ob -j^xj dMn(M,z)p(M ob \\M), (3.1) 

2 1 l,m ^ 

where dV/dz is the comoving volume element per unit redshift and solid angle, and n(M, z) 
is the halo mass function. In this notation M° b _ = M m j n (z) represents the minimum value 
of the observed mass for a cluster to be included in the survey, and it is determined by 
the survey selection function (see § 3.4). The integral over the observed mass is computed 
within bins having width AlogM = 0.2, extending from M m \ n (z) to 10 15,8 /i _1 Mq. For 
n(M,z) we adopt the expression provided by [52], with mass function parameters obtained 
for overdensity A = 200 with respect to the mean density of the universe (see their Table 
2). Moreover, to take into account massive neutrino effects, we follow the prescription used 
by many authors [e.g. 53-55] neglecting the weakly clustering neutrino component when 
calculating the halo mass (M ha io = 47rr 3 p/3, with p 

— Pm ~ Pv)- Many other calibrations 
of the halo mass function from simulations have been presented by several authors [e.g. 56- 
59]. However, for the purposes of this work the choice of the best-calibrated mass function 
has a minor impact. Indeed the forecast errors depend primarily on the number of cluster 
expected for a given cosmological model, which is far more sensitive to the exponential shape 
of the mass function rather than to the calibration details of the mass function. The factor 
p(M ob \\M) takes into account the uncertainties that a scaling relation introduces in the 
knowledge of the true cluster mass. Following the prescription of [60], p(M ob \\M) gives the 
probability of assigning to a cluster of true mass M an observed mass M ob , as inferred from 
a giving scaling relation. Under the assumption of a lognormal-distributed intrinsic scatter 
around the nominal scaling relation with variance crf nM , the probability can be written as: 

p(M° b \\M) = -P[-* 2 (^)] , (3 .2) 

V ( 27T(J Im) 
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where 

, nh , lnM ob - B M -\nM , . 

x(M ob ) = . (3.3) 

Here the parameter Bm represents the fractional value of the systematic bias in the mass 
estimate. Moreover, according to [61] we assume the following parameterization for red- 
shift dependencies of the halo mass bias and variance (we do not consider a possible mass 
dependence of these parameters): 

B M (z) = £ M ,o(l + z) a 

0"lnM(» = 0"lnM,o(l + zf . (3.4) 

In our formalism, we have four nuisance parameters, -£>M,(b %moi a an d ft, which can be 
allowed to vary along with the other cosmological parameters during the forecast procedure 
(see § 4.4). 

Including Eq. 3.2 into Eq. 3.1 it follows that: 

Af2 f Zl+1 dV f°° 
Ni, m = — J dz ~ d JJ dMn ( M > z ) x [erfc(x lim ) - erfc(xi, m+1 )] , (3.5) 

where erfc(x) is the complementary error function. In the following, we assume that the errors 
on cluster redshift measurements can be neglected. For Euclid data, the photometric redshift 
measurements will be calibrated using a combination of the spectroscopic survey and ground- 
based visual bands photometry, with an expected limiting precision of a(z) ~ 0.05(1 + z) [62]. 
While this error refers to the photometric redshift of a single galaxy, the redshift of a cluster 
identified in the photometric survey should be reduced by a factor iV 1 / 2 , where N is the 
number of galaxies assigned to the cluster. 

3.2 Cluster power spectrum 

In order to include information from the clustering of galaxy clusters, we calculate the aver- 
aged cluster power spectrum within a given redshift interval using the expression: 

P?(k,z)= Jz < * \ (3.6) 



where n(z) = J °° dMn(M, z) x erfc(xi iin= o) is the comoving number density of clusters that 
are included in the survey at redshift z [e.g. 63]. The cluster power spectrum P cl (k,z) is 
expressed in terms of the underlying matter power spectrum P(k, z) according to P cl (k, z) = 
bcg(z) P(k, z); the term of proportionality b e ft is the cluster mass function averaged linear 
bias, defined as: 

J °° dMn(M, z) erfc(xi, m=0 ) b(M, z) 
beS[Z) J^dMn(M,z)evfc(x hm=0 ) ' 

For the linear bias of dark matter halos, b(M, z), we adopt the fitting function of [64] for over- 
density A = 200 (see their table 2). The linear matter power spectrum P(k,z) is computed 
with the publicly available software package CAMB [65] , which takes correctly into account the 
effect of massive neutrinos also in a mild non- linear regime [66]. In order to avoid contribu- 
tion to the matter power spectrum and scale-dependent bias introduced by non-linearities, we 
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Figure 2: Wavenumber dependence of the relative error of the cluster power spectrum, 
<T pel /P cl defined as in Eq. 3.10, at four different redshift: z = 0.2,0.6, 1.0, 1.6, from bottom 
to top curves, respectively. 

do not include in our analysis modes with wavenumbers larger than fc max = 0.1 Mpc -1 [67]. 
Although massive free-streaming neutrinos mainly affects the power spectrum at small scales 
(see section 2), a value of A; max larger than ~ 0.3 Mpc would not increase significantly the 
sensitivity of the survey. Indeed, given the level of Poisson noise associated to the cluster 
distribution (see Eq. 3.10, and Fig. 2), higher frequency modes are not adequately sampled 
and, therefore, adding them to the analysis does not add significant information [61]. As for 
the minimum value of the wavenumber we impose k m i n = 0.003 Mpc; again, using a smaller 
value of k m [ n does not change the final results, since extremely large scales are not sampled 
by the surveys. Finally, in our analysis we neglect the correction to the power spectrum due 
to redshift space distortion effects. 

3.3 Forecasting 

The forecast is based on the Bayesian inference technique, for which a likelihood function of 
the mock data is first constructed and then sampled in order to estimate the marginalized 
probability distribution of the parameters. To explore the parameter space by means of 
Monte Carlo Markov Chains we use the publicly available code CosmoMC 8 [68], where we 
included a module for the calculation of the cluster number counts and power spectrum 
likelihoods. 

Our most general parameter space is: 



where the first six parameter, which define the standard ACDM model, are: the physical 
baryon Wb = ^b^- 2 an d c °ld dark matter uj c = Q, c h 2 densities, the ratio (multiplied by 
100) between the sound horizon ant the angular diameter distance at decoupling S , the 

8 http : / /cosmologist . info/cosmomc / 



= (w b ,W c ,e s ,T,n s ,log[10 A s ], f^Nefi^^W, B M ,0, <JlnM,0, a, P) 



(3.8) 
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Table 1: Fiducial parameter values. 
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3.046 





-1 





0.45 









reionization optical depth r, the scalar spectral index n s and the amplitude of initial power 
spectrum ^4 S . Besides these parameters we performed several forecasts for different extensions 
of the minimal cosmological model, by fitting (along with the other parameters): the neutrino 
density fraction /„ = £l u /£l c , the effective number of neutrino species N c q, the spatial 
curvature f2k and the dark energy equation of state parameter w. Finally, in order to assess 
the effect of the uncertain knowledge of the mass-observable relation we consider also the 
case in which the four nuisance parameter are treated as fitting parameters to be determined 
along with the cosmological ones. 

Throughout this paper, our reference model is chosen to be a flat ACDM model with 
three neutrino species. The fiducial ACDM parameter values are listed in Table 1, consistently 
with the WMAP-7+BAO+i?o best-fit model by [69]. These fiducial parameter values are 
also consistent with the latest WMAP 9-yr best- fit model [9]. 

As for the nuisance parameters, clusters mass within Euclid survey will be estimated 
using photometric richness as a cluster mass proxy. Since the scaling relation between richness 
and mass will be calibrated using weak lensing measurements, which directly map the matter 
distribution of the cluster probing its gravitational tidal field, we consider £?m,o — as 
a reference value for the mass bias. For the intrinsic scatter we assume <Xi n M0 = 0-45 
as estimated by [39] by demanding consistency between available weak lensing and X-ray 
measurements of the maxBCG clusters, and the X-ray luminosity-mass relation inferred from 
the 400d X-ray cluster survey. The intrinsic scatter has the effect of increasing the number of 
clusters included in the survey. Indeed, the number of low-mass clusters that are up-scattered 
above the survey mass limit is always larger than the number of rarer high-mass clusters which 
are down-scattered below the same mass limit [e.g. 70]. Because so far there are no evidences 
for the evolution of the nuisance parameters we adopt a = and j3 = as reference values, 
thus making the minimal assumption of constant bias and scatter with redshift. When the 
nuisance parameters are left free we consider two cases: one with strong prior on the evolution 
parameters, with a and f3 not allowed to vary with respect to their reference value, and the 
other one with no prior knowledge of their value. The latter turns out to be a conservative 
assumption in view of the large number of clusters for which mass measurements from weak 
lensing will be available from Euclid [62]. As such, the corresponding uncertainties expected 
on cosmological parameters should be regarded as upper limits of the error introduced by the 
uncertainties in the relation between cluster richness and mass. The cluster power spectrum 
and number counts of the mock data are assumed to be equal to the theoretical cluster power 
spectrum and number counts of the fiducial model. 

Since we are interested only in parameter error estimation, we define our likelihood 



- 8 - 



functions C of the observable O as 



xls= -21n£ 



E — 



obs 



(3.9) 



in such a way that xts is equal to zero for the fiducial parameter values. In the previous 
equation 0°^ s denotes the observed cluster power spectrum, P°^ s (ki, Zj) (number counts, 
iV obs (Mj, Zj)), while O^j 1 is the theoretical cluster power spectrum (number counts) of Eq. 3.6 
(Eq. 3.5). The statistical error associated to the observed galaxy cluster power spectrum in 
a bin centered on (ki,Zj) [71] reads: 



2 

O pel 

ij 



(2vr)^(P c r(^j)) 
V 8m (zj)kfAk 



I 2 



1 + 



obs / 



(3.10) 



where V SUT (zj) is the comoving survey volume within the redshift bin centred on Zj, and Ak 
is the size of the bins in wavenumber space. In this way, constraints at redshift z are mostly 
contributed by wave-numbers k, which maximize the product n(z)P°]° s (k, z) (see Fig. 2). 
The average cluster power spectrum is computed by integrating over redshift intervals having 
constant width Az = 0.2. This choice of binning represents a compromise between the need 
of extracting the maximum amount of information from clustering evolution and request of 
negligible covariance between adjacent z-intervals [e.g. 72]. Indeed, the definition of Eq. 3.9 
holds only if the contribution from different redshift slices carry statistically independent 
information. As for the statistical error of the observed number counts for a given mass and 
redshift bin centered in (Mj, zj), we consider only the Poissonian noise, a^ cl = N ohs (Mi, zj), 

neglecting the contribution from sample covariance (i.e. cosmic variance), which accounts 
for the clustering of clusters due to large scale structure (for a more rigorous definition of the 
number counts covariance matrix see e.g. [73]). 

Because the full parameter space is quite large and some parameters are poorly con- 
strained by LSS observations, we perform our forecast combining the Euclid-like cluster 
catalog with Planck-like data. The mock CMB TT, EE and TE power spectra have been 
simulated following the procedure of [41] according the specifications presented in the Planck 
Blue-Book [74, page 4, Table 1.1] based on 14 months of observations, using the three fre- 
quency channels with the lowest foreground levels at 100GHz, 143GHz and 217GHz, and 
a sky fraction of / s k y = 0.80. In order to avoid problems with foreground signal, beam 
uncertainties, etc., we cut-off the spectra, set l m ax — 

2000. 



3.4 Characteristics of the survey 

Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, planned 
for lunch in 2019. Thanks to its three imaging and spectroscopic instruments working in the 
visible and near-infrared bands, Euclid will cover 15,000 square degrees of extragalactic sky 
with the wide survey, thereby providing high-quality images for from visual imaging for more 
than a billion galaxies, accurate photometric redshifts from near-IR imaging photometry 
(in combination with ground-based data) for about 2 x 10 8 galaxies and about 5 x 10 7 
spectroscopic redshift at z > 0.7 from near-IR slitless spectroscopy. 

The most efficient method to build the Euclid galaxy cluster catalog relies on the anal- 
ysis of the photometric data. To predict cosmological constraints from the expected sample 
of galaxy clusters, we use the analytic selection function adopted in the Euclid Red-Book [62] 
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Figure 3: The cumulative cluster redshift distribution as predicted by the reference cosmo- 
logical model and the reference values for the mass nuisance parameters (see Table 1), for 
the Euclid cluster survey. 

to forecast the contribution of the cluster survey to the cosmological constraints. The com- 
putation of the selection function is based on using the luminosity function of cluster galaxies 
to compute the number of galaxies expected within i?200e 9 down to the Hab = 24 magni- 
tude limit reached in the photometric survey, as a function of the cluster mass. Comparing 
this number with the expected fore-/back-ground galaxies, and including the dilution effect 
associated to photometric redshift errors, one can compute the expected S/N with which a 
cluster of a given mass, obtained from a calibrated mass-richness relation, can be detected 
as a function of redshift. Assuming S/N = 3 for the fiducial limiting signal-to-noise for a 
reliable cluster detection turns into a selection function which provides Mii mi 2ooc(z) defined 
as the limiting mass within i?200e for a cluster to be included in the survey. As a result, one 
finds Mii mj 2ooc(-z) — 1.6 x 1O 14 M at z > 0.5, while decreasing at lower redshift, reaching 
~ 5 x 1O 13 M at z = 0.2. Finally, to compute Eq. 3.5 we need to convert Mi im)2 ooc( 2; ) to 
-^Hm,200m(^) ~ the limiting mass within a radius encompassing an overdensity equal to 200 
times the mean density of the Universe - consistently with the chosen halo mass function 
(see § 3.1). To this end we fallow the recipe given in [73], assuming a NFW profile [75] as 
halo density profile and using their fitting formula (Cll). 

For the sky coverage we adopt the required (covered) area for the wide survey AO = 
15,000deg 2 , while the cluster number counts and power spectrum are evaluated in redshift 
bins of width Az = 0.2 between z = 0.2 and z = 2. From Eq. 3.5, using this survey specifics 
and fiducial parameter values, we expect that Euclid will find of order 1.5 x 10 5 cluster with 
a S/N better than 3, between z = 0.2 and z = 2.0, with ~ 4 x 10 4 having z > 1 (see Fig. 3). 



9 Here i?200c is defined as the radius encompassing an average density equal to 200 times the cosmic critical 
density at a given redshift 
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Table 2: Constraints on E m v for ACDM+m„ model from Planck, cluster power spectrum 
(P cl — only), Euclid-Cl (cluster number counts and power spectrum) data, and the combina- 
tion of the two data sets Euclid-Cl+Planck. Because the parameter r is not constrained by 
Euclid data, when CMB measurements are not included r is kept fixed to its fiducial value. 



Model 






ACDM+m„ 




Data 


Planck 


P cl — only 


Euclid-Cl 


Euclid-Cl+Planck 


E™. [eV] 


68% CL < 0.41 
95% CL < 0.74 


< 0.60 

< 1.37 


< 0.17 

< 0.35 


< 0.017 

< 0.031 



4 Analysis and results 

Having defined the reference cosmological model and the specifics of Euclid survey we now 
present forecast errors on neutrino parameters for various extensions of the minimal ACDM 
model. For each of the cases that we describe here below, we run four independent chains, 
requiring the fulfillment of the Gelman & Rubin [76] criteria with R— 1 ^ 0.03 as convergence 
test. 

4.1 Three massive neutrinos: ACDM+m^ 

We start by considering the scenario with three degenerate massive neutrino species. In 
Table 2 we report the 68% and 95% CL bounds on E m v derived from different data sets: 
Planck only, cluster power spectrum only, cluster number counts and power spectrum (here- 
after Euclid-Cl), and the combination of Planck and Euclid-Cl. When only cluster power 
spectrum data are considered we obtain a quite loose 2a upper limit on E m f °f l-37eV. 
Otherwise, the information contained in cluster number counts alone is unable to constrain 
the total neutrino mass, but it greatly improves the error on E m v once added to cluster 
power spectrum data, mainly thanks to the thigh constraints provided on ag (see left panel 
of Fig. 4). Specifically, the upper limit for E m u shrinks by a factor ~ 4 to 0.35 eV (95%CL). 
This error is comparable to the present constraints obtained combining CMB and LSS probes, 
and of the same order of magnitude of the error expected for Planck. Regarding parameter 
degeneracies for the galaxy cluster dataset, the total neutrino mass is correlated with all the 
cosmological parameter affecting the galaxy power spectrum shape (i.e. O m , erg, n s ; see red 
contours in Fig. 4 and Fig. 5). 

The main power of constraints in cosmological parameters indeed originate from the joint 
analysis of galaxy cluster and CMB datasets. In this case the error on E m v i s reduced to 
31 meV; an improvement of more than one order of magnitude that would allow a 2a detection 
of the total neutrino mass even in the minimal normal hierarchy scenario (E m u — 0.05 eV). 
The reason for such an improvement can be easily understood by looking at the right panel 
of Fig. 4 which shows the 68% and 95% confidence regions in the i^2m u — as) plane from 
Planck data, Euclid-Cl data and the combination of the two. Taken independently, the 
CMB and galaxy cluster data exhibit significant degeneracies in this plane, but the nearly 
orthogonal degeneracy directions allow their combination to provide tight constraint on these 
parameters, and in particular on the neutrino mass. When Planck priors are added to the 
Euclid-Cl constraints, all degeneracies are either resolved or largely reduced (see blue contours 
in Fig. 4 and Fig. 5). Similar levels of sensibility on E m u are also expected combining Euclid 
galaxy or cosmic shear power spectrum measurements with Planck CMB data [e.g. 77-79] 
and from the combination of Planck SZ cluster survey and Planck CMB data [80]. 
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Figure 4: The 68% and 95% CL contours in the i^m^-ag) plane for a ACDM+m„ model. 
Left panel: contours from cluster power spectrum (large contours; P c '-only) and the combi- 
nation of cluster power spectrum and number counts (small contours; Euclid-Cl). The insert 
plot shows a zoom of the confidence contours given by the Euclid-Cl dataset compared with 
the contours obtained from the Fisher Matrix technique using the same dataset. Right panel: 
contours from Planck (green), Euclid-Cl (red) and Planck -l-Euclid-Cl (blue) datasets. The 
insert plot shows a zoom of the confidence contours obtained from the Planck -l-Euclid-Cl 
datasets. 
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Figure 5: The 68% and 95% CL contours in the — (il m ,n s ) planes for a ACDM+mi/ 

model, from Euclid-Cl (red contours) and the Planck -l-Euclid-Cl (blue contours) datasets. 
When only Euclid-Cl dataset is used the parameter r, which is not constrained by this data, 
is kept fixed to its fiducial value 0.085. 



4.2 Varing jV eff : ACDM+m v +N eH 

We now explore the scenario with massive neutrinos and N e g effective number of neutrino 
species. Again, we distribute the sum of neutrino masses equally among three active species 
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Table 3: Constraints on ^ m„ and iV e g for ACDM+m^+A^g model. 



Model 




ACDM+m^+iVcff 


Data 




Planck Euclid-Cl+Planck 


m u [eV] 


68% CL 
95% CL 


< 0.42 < 0.022 

< 0.78 < 0.040 


iVcff 


95% CL 


< 3.36 < 3.14 



{N v = 3), and we treat additional contribution to iV e g as massless, such that N e g = 3 + AiV c g, 
with the prior AA^ e jj ^ 0. While the choice of keeping N v fixed does not affect constraints 
from CMB measurements (what matters is N e s ) , it could change the sensitivity to m v 
and A^ e fj based on galaxy clusters data. Indeed, changing N v would change the mass of each 
massive neutrino and thus its free-streaming lenght beyond which the power spectrum is 
suppressed (see Eq. 2.3). We checked the case with fixed AN e g and N v as free parameter 
and we find no qualitative changes in our results. This means that the data are not sensitive 
to the exact position of the break in the power spectrum induced by free-streaming neutrinos, 
and thus to the neutrino mass hierarchy. Table 3 shows the joint constraints on the sum of 
neutrino masses and on the effective number of neutrino species from Planck data alone 
and the combination of Planck and Euclid-Cl datasets. Looking at Planck data alone, the 
quality of the constraints on ^ m u are nearly unchanged from the single-parameter extensions 
discussed earlier, as it would be expected for independent parameters. Indeed, Yl m v an d AT e ff 
are constrained by different features in the CMB spectra: the early integrated Sachs- Wolfe 
effect for Y rrL v, the damping scale and the position of the acoustic peaks for iV e g [see e.g. 
9, and references therein ]. However, since both Yl m f an d -^eff have similar effects on the 
matter power spectrum (see § 2), the correlation of the two degrades the upper bound on the 
sum of neutrino masses inferred from Euclid-Cl+Planck data by ~ 30% to ^ m u < 0.040[eV] 
at 95%CL. With this accuracy, it would still be possible a 2 a detection of neutrino masses in 
the minimal normal hierarchy scenario. Constraining N e g is mainly achieved through CMB 
measurements of the redshift of the matter-radiation equality z eq and the baryon density 
fl^h 2 . However, keeping z eq and Q.\ i h 2 fixed as N e g increases can be achieved by increasing 
the cold dark matter density £l c h 2 , which displays a large correlation with N e $ [81]. Euclid-Cl 
data alone are unable to provide constraints on the number of effective species. However, the 
inclusion of clusters dataset allows to significantly improve the measurements of £l m h 2 (by 
constraining as), thus reducing the 2a error on the effective number of neutrino by a factor 
larger than 2.5 from 0.36 to 0.14 (see right panel of Fig. 6). After the inclusion of Euclid-Cl 
data N c ff still exhibits strong degeneracies with many cosmological parameters (e.g. VL m h 2 , 
Hq and n s ), and a correlation of ~ 0.5 with ^ m u (see left panel of Fig. 6). 

4.3 Extended models: zuCDM+m^+iVeff and curved Universe 

Next, we consider how the constraints on ^ m v and N e f[ are affected when additional degrees 
of freedom are introduced in the cosmological model. The effect on (Y m u, N e s) of adding 
these degree of freedom to the ACDM+m^+A^eg model are shown in Fig. 7 and listed in 
Table 4. We start by considering a constant dark energy equation of state (w ^ —1)- Neutrino 
properties (Y rn u , N e g) and w are generally degenerate because they can both affect the shape 
of the matter and CMB power spectra [e.g. 82]. Looking at Fig. 8 (a), we indeed see this 
degeneracy in the plane N e ^—w, which displays a correlation of ~ 0.5, whereas the parameters 
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Figure 6: The 68% and 95% CL contours in the iV e g — (Q m ,E m ^) planes for a 
ACDM+m^+iVeff model, from: Planck (green contours) and Planck+Euclid-Cl (blue con- 
tours) data. 

Table 4: Constraints on E m v an d N e s for the two parameter extensions (w, Cl^) from 
Euclid-Cl+Planck datasets. 



Data 




Planck+Euclid-Cl 


Model 




wCBM+m u +N cG ACDM+m^+iVeff+Ok 


[eV] 


68% CL 


< 0.024 < 0.024 


95% CL 


< 0.046 < 0.046 


iVeff 


95% CL 


< 3.16 < 3.17 



w and E m v show almost no correlation. The Euclid clusters catalog, probing the evolution 
of the LSS up to z ~ 2, will be able to put tight constraints on the dark energy equation 
of state; we find for the combination of Planck and Euclid-Cl data: —1.011 < w < —0.987 
(95%CL). Given the small uncertainty on w the constraints on neutrino mass and effective 
number of species are only slightly degraded when w is allowed to vary; the 95% CL upper 
limit for iV e fj is relaxed from 3.14 to 3.16 due to the degeneracy with w. Whereas, the 95% 
CL upper limit for undergoes only a small degradation from 0.040 eV to 0.046 eV, 

caused by the weak constraints on parameters that are correlated with E m v induced by 
the extension of the parameter space. Secondly, we relax the prior on the curvature of the 
universe by considering the case ACDM+m I/ +A' e ff + $7^. Since current data do not support 
departures from the flat ACDM model either through ^ or w ^ — 1, we introduce these 
parameters separately. From the combination of Planck and Euclid-Cl datasets we obtain, 
for the curvature parameter, the following constraint: —0.0024 < O^ < 0.0024 (95%CL). 
The spatial curvature mainly affects the expansion rate via the Friedmann equation, as well 
as the total neutrino mass and number of effective species do. As it can be seen in Fig. 8 
(b), this results in a correlation with both Y2 m v anc l -^cff of the order of ~ 0.5 and ~ 0.6, 
respectively. Despite these quite large degeneracies with i\, the small error associated to the 
curvature parameter leads to a slight relaxation of the constraints on neutrino properties: 
the upper limit for neutrino mass degrades by ~ 10%, passing from 0.040 eV to 0.046 eV 
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Figure 7: The marginalized one-dimensional posteriors for ^ m v (left) and iV e fT (right) for 
different parameter extensions from the combination of Euclid-Cl and Planck datasets. 

Table 5: Constraints on ^ m v an d N^r f° r ACDM+m^-)- iV c fj models with free nuisance 
parameters. 



Data 




Planck+Euclid-Cl 


Model 




ACDM+m^+Afcff+nuis ACDM+m^+A^ff+nuiscv 


E m v [eV] 


38% CL 


< 0.049 < 0.031 


35% CL 


< 0.083 < 0.056 




35% CL 


< 3.18 < 3.16 



(95%CL), while the 2a error on N cS shift from 0.14 to 0.17, a 20% degradation. 

Thus, in both cases, the parameter extension entails a relaxation of the constraints on 
m v and N e g; nonetheless, given the high accuracy with which w and are expected to 
be measured, the survey would still allow a 2a detection of neutrino mass in the minimal 
normal hierarchy scenario and reveal the presence of possibles extra relativistic species. 

4.4 Nuisance parameters 

Finally, to assess the effect of an uncertain knowledge of cluster masses on ^ m„ and N e g 
constraints, we treat the ACDM+m^ + N e Q case with the four nuisance parameters as fitting 
parameters, following the so-called self-calibration method [e.g. 51, 60, 63]. The results are 
listed in Table 5. We start with the the over-conservative assumption of no priors on all 
the four nuisance parameters (ACDM+m^+A r c fj-|-nuis model). The uncertainties on scaling 
relation parameters compromise our ability to recover the halo mass function from cluster 
data, thus reducing the cosmological information achievable from cluster number counts. 
This results in a larger error for the parameters that are primarily constrained by cluster 
number counts, in particular for as, the normalization of the power spectrum. Looking 
at Fig. 9 (left panel), the constraints on ag are relaxed by a factor of ~ 10 compared to 
the ACDM+m^+A'eff model, and the parameter recovers a large degeneration with m v 
of the order of ~ 0.65. This effect, along with the degradation of other parameters errors 
(e.g. a(Q m )), entails a relaxation of the upper limit for by a factor larger than two, 

from 0.040 eV to 0.083 eV. With these loose constraints, in the case of minimal normal 
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Figure 8: Joint two dimensional marginalized constraints on w (a) and fik (b) against 
(^2m Ul N c g) at 68% and 95% CL. The confidence regions are respectively for the extended 
parameter space wCDM+m u +N c s and ACDM+m 1/ +-W e flE + fik, using data from Euclid- 
Cl+Planck. 



hierarchy scenario, it would not be possible to have a two a detection of neutrino mass. 
Because the constraints on neutrino mass from cluster number counts relay on the evolution 
of the high-mass end of the mass function, ^ m u is rather degenerate with a and /3, the two 
nuisance parameters which control the evolution of the systematic bias and intrinsic scatter 
(see Eq. 3.4). To emphasize the role played by the uncertain redshift evolution of the nuisance 
parameter on the determination of ^ m„ we show in Fig. 9 the contours for a model with 
a and /3 kept fixed (ACDM+m 1 ,+.ZV e flH-nuis e .u model). In this case the degradation of the 
total neutrino mass constraints with respect the ACDM+m^+A^g model is only of ~ 40%, 
from 0.040 eV to 0.056 eV. In other words, an accurate knowledge of the redshift evolution 
of the nuisance parameter improves the 2a upper limit of m v by ~ 33% compared to the 
previous case with no prior on the nuisance parameters. 

Likewise, the forecast error on iV e fj is influenced by the loss of constraining power of the 
cluster number counts data, even if to a lesser extent than the bounds on Y2 m vi since the 
constraints on N e g are primarily contributed by CMB measurements. The 2a upper limit 
shifts from N c g < 3.14 to iV e g < 3.16 and N e g < 3.18, in the model with strong evolution 
prior and free nuisance parameters, respectively. In these case the degradation is mainly due 
to the larger error associated to Cl m h 2 , which is highly degenerate with N e g as explained 
in section 4.2 and shown in Fig. 9 (right panel). We remind that the results for the model 
with no prior have to be regarded as an upper limit on the error introduced by the uncertain 
knowledge of the scaling relation; nevertheless, these results highlight the importance of 
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Figure 9: Joint two dimensional marginalized constraints on the planes (Y^ m u — fs) an d 
(N cS - n m h 2 ) at 68% and 95% CL from Euclid-Cl+Planck data. The confidence regions are 
for the ACDM+m I/ +iV e ff model discussed in 4.2 (small blue contours) and the two extended 
model with nuisance parameters: all nuisance float (larger light violet contours) and fixed 
evolution parameters a and (3 (dark violet contours). 

having robust calibration of the scaling relation, and in particular of their evolution with 
redshift, to fully exploit the cosmological information contained in the Euclid cluster catalog. 

5 Conclusions 

In this paper, we presented forecasts on the capability of a future photometric galaxy clus- 
ter survey such as Euclid, in combination with Planck-like data, to provide constraints on 
neutrino properties. Specifically, we rely on two observables: the cluster number counts and 
their power spectrum. Our analysis is based on the Markov Chain Monte Carlo methods 
rather than the Fisher Matrix technique, which results in more reliable error bars. We start 
by considering a reference ACDM model in agreement with the recent results of WMAP 9-yr. 

In order to study possible degeneracies with ^2m u , besides the ACDM model with 
massive neutrino, we also consider models with iV e ff effective number of relativistic species, 
a constant dark energy equation of state w and curvature. Following the self-calibration 
approach, along with the other cosmological parameter, we decide to explore also the effect 
of leaving free the nuisance parameters that describe the relation between cluster optical 
richness and mass, its scatter and redshift evolution. 

Our results can be summarized as follows: 

• From the combination of Euclid number counts and clustering data we obtain a 2a 
upper limits for the total neutrino mass of ^2m u < 0.35 eV, comparable with present 
constraints from the combination of CMB and LSS probes [e.g. 6, 9, 37]. When Planck 
data are added to the Euclid-Cl ones the error on Yl m u is reduced by a factor larger 
than 10 to m v < 0.031 eV. With this accuracy the total neutrino mass could be 
detected at 2a level even in the minimal normal hierarchy scenario. The large improve- 
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ment is due to the different degeneracies present between Euclid and Planck that are 
broken once the two experiments are combined. 

• Because the effective number of neutrino spices is degenerate with the sum of neu- 
trino masses, varying iV e g entails a relaxation of 2a error bars on 52 m u by ~ 30% in 
the Planck+Euclid-Cl case. Still, the 2a error is lower than the minimum neutrino 
mass admitted by neutrino oscillation experiments. The Euclid-Cl dataset is unable 
to constraints N e g by itself, but improves the 2a upper limits on N e g from 3.36 using 
Planck-only, to 3.14 in the Planck+Euclid-Cl case. The improvement is mainly due to 
the tighter constraints on VL m h 2 provided by the Euclid-Cl datasets. 

• In models with varying w or £1^ the 2a error on 52 m u is relaxed only by ~ 10%. In both 
cases the high accuracy with which w or are constrained by the Planck+Euclid- 
Cl data prevents the error on 52 m u from being largely degraded. As for N e s, the 
parameter shows a correlation of the order of ~ 0.5 with both w and Ok; which shifts 
the 2a upper limit for iV c fj to 3.16 and 3.17, respectively. 

• When nuisance parameters are considered in a conservative way (no prior) our ability 
to recover the halo mass function from cluster data is compromised. The degradation 
of cosmological information results in a ~ 2 times larger 2a error for neutrino masses 
(52 m v < 0.083 eV) and a degradation of ~ 30% of the 2a error on the effective num- 
ber of neutrinos (N e g < 3.18). In this case the accuracy would not be sufficient for 
detecting the total neutrino mass with good significance in the minimal normal hierar- 
chy scenario. Whereas, assuming a perfect knowledge of the redshift evolution of the 
nuisance parameters we partially recover the informations contained in cluster number 
counts data. In this case the 2a upper limit for 52 m v i s degraded only by 40% to 
52 m v < 0.056 eV, while the 2a error on the effective number of neutrinos degrades by 
~ 15% to N cS < 3.16. 

It is worth reminding that in our analysis we did not include the effect of redshift space 
distortions in the distribution of galaxy clusters induced by peculiar velocities. This effect 
should be in principle included when forecasting the cosmological constraining power of future 
cluster surveys; indeed, as demonstrated by [51], the inclusion of redshift space distortions 
carries significant cosmological information through the growth rate of density perturbations. 

As a concluding remark, we emphasize once again the importance to provide an accurate 
calibration of the scaling relation between the observable quantity on which cluster selection 
is based, optical richness in this case, and cluster mass. Thanks to the exquisite imaging 
quality expected for the Euclid survey, weak lensing masses for individual objects will be 
available for a significant fraction of the clusters identified from the photometric selection. 
At the same time, staking analysis will provide an accurate calibration of the relation between 
weak lensing masses and richness. For this reason, the Euclid cluster survey will represent a 
powerful complement to galaxy clustering and cosmic shear analyses to constrain cosmology 
through the growth of perturbations. 
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